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ABSTRACT 

The Sc galaxy M 99 in the Virgo cluster has been strongly affected by tidal interactions 
and recent close encounters, responsible for an asymmetric spiral pattern and a high 
star formation rate. Our XMM-Newton study shows that the inner disk is dominated 
by hot plasma at kT « 0.30 keV, with a total X-ray luminosity w 10 41 erg s _1 in 
the 0.3-12 keV band. At the outskirts of the galaxy, away from the main star-forming 
regions, there is an ultraluminous X-ray source (ULX) with an X-ray luminosity « 
2 x 10 40 erg s _1 and a hard spectrum well fitted by a power law of photon index 
r w 1.7. This source is close to the location where a massive Hi cloud appears to be 
falling onto the M99 disk at a relative speed > 100 km s _1 . We suggest that there 
may be a direct physical link between fast cloud collisions and the formation of bright 
ULXs, which may be powered by accreting black holes with masses ~ 100M©. External 
collisions may trigger large-scale dynamical collapses of protoclusters, leading to the 
formation of very massive (> 200 Mq) stellar progenitors; we argue that such stars 
may later collapse into massive black holes if their metal abundance is sufficiently low. 
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1 INTRODUCTION 

The Sc galaxy M 99 (NGC 4254), located a distance of about 
17 Mpc (Tully 1988), is the brightest spiral in the Virgo 
Cluster (Mb = —20.8), with a gas mass « 5 x 1O 9 M0 and 
a kinematic mass « 10 Mq (Vollmer, Huchtmeier & van 
Driel 2005; Phookun, Vogel & Mundi 1993). It shows a pe- 
culiar spiral structure 1 , with one arm less tightly wound 
and much brighter than the other two. Another unusual 
feature of M 99 is the presence of a string of H I "blobs" 
to the south and west of the stellar disk, and a low-surface- 
density HI gas tail to the north-west; the total gas mass 
in those extra-disc structures is ~ 2 x 10 s Mq (Phookun et 
al. 1993). It is likely that there is a physical connection be- 
tween the unusual spiral structure and the disturbed gas 
distribution. In one scenario (Vollmer et al. 2005), the spi- 
ral structure was affected by a close encounter with another 
Virgo Cluster galaxy « 280 Myr ago; the surrounding HI 
clouds are due to ongoing, face-on ram-pressure stripping, 
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1 As a matter of historical curiosity, it was the second galaxy 
ever in which a spiral pattern was identified (Rosse 1850), with 
the "Leviathan of Parsonstown" , a few months after a similar 
discovery for M51. 



and are mostly moving away from the galaxy. An alterna- 
tive scenario (Phookun et al. 1993) suggests that the HI 
clouds and tail are tidal debris of an entity that was at least 
partly disrupted in a close encounter with M99; that debris 
is now infalling and merging with the galactic disk. Simu- 
lations show (Bekki, Koribalski & Kilborn 2005) that tidal 
debris could also have been stripped from the outer disk 
of M 99 itself, by the Virgo Cluster potential as the galaxy 
crossed the central region of the cluster; part of the stripped 
gas would later fall back onto the same galaxy. In this sce- 
nario, heavy gas infall onto the M 99 disk would be respon- 
sible for the lopsided spiral structure (Phookun et al. 1993; 
Bournaud et al. 2005). 

Intriguingly, a large H I cloud ( VIRGOHI 21) , with a gas 
mass ss 2 x 10 s M© but without any associated stars, located 
~ 120 kpc to the north-west of M99, has recently been dis- 
covered (Davies et al. 2004). One interpretation (Minchin et 
al. 2005, 2006) is that it is an old, bound "dark galaxy" that 
has never formed stars because of its low hydrogen surface 
density. In this scenario, VIRGOHI 21 would be dark-matter 
dominated, with a kinematic mass ss 1O 11 M0; if so, it would 
likely be the same galaxy responsible for the suspected close 
encounter with M 99 and possible tidal gas stripping. Alter- 
natively, VIRGOHI 21 could itself be simply another, larger 
piece of tidal debris (without a dark matter halo), stripped 
from the outer H I disk of M 99 by the Virgo Cluster poten- 
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tial or during a close encounter with another galaxy (Bekki 
et al. 2005) 

Gas-rich galaxies with recent tidal interactions gener- 
ally show particularly active star formation, and M 99 is no 
exception, with a rate ~ lOM© yr _1 inferred from its Ha lu- 
minosity (Kennicutt et al. 2003), or m 1O-2OM0 yr _1 from 
an optical/near-IR photometric study (Gonzales & Graham 
1996). As a comparison, this is twice as high as in M 82 (Ken- 
nicutt et al. 2003), and three times the current total rate in 
our Local Group (Hopkins, Irwin & Connolly 2001). An- 
other phenomenon often associated with tidal interactions 
and high star formation is the presence of ultraluminous X- 
ray sources (ULXs) with apparent isotropic X-ray luminosi- 
ties ~ 10 39 -10 40 erg s _1 (for recent reviews, see King 2006; 
Fabbiano & White 2006; Colbert & Miller 2004). They are 
interpreted as accreting binary systems, an order of mag- 
nitude more luminous than the Eddington limit for typi- 
cal stellar-mass black holes (BHs) in the Local Group. Two 
ULXs were detected by ROSAT/YiRl within the D25 region 
of M 99, on 1997 June 30 (Colbert & Ptak 2002; Liu & Breg- 
man 2005). In this paper, we use XMM-Newton to have a 
better understanding of the X-ray properties of M99, and 
in particular of its brightest ULX. We compare radio and 
X-ray information to investigate the possible connection be- 
tween ULX formation, star formation and collisional events 
in this galaxy. Finally, we outline a possible general scenario 
for the nature of ULXs, to be tested by future studies. 



2 XMM-NEWTON STUDY 
2.1 Data analysis 

M 99 was observed with all instruments on-board XMM- 
Newton on 2003 June 26 (Revolution 651; observation ID 
0147610101); the thin filter, full-frame mode was used for 
the European Photon Imaging Camera (EPIC) detectors. 
We processed the Observation Data Files with standard 
tasks in the Science Analysis System (SAS), version 6.5.0. 
Unfortunately, most of the X-ray observation was disrupted 
by strong solar flares: after inspecting the background 
lightcurves, we retained a good live-time interval of 17.0 
ks for the pn and 20.5 ks for the MOSs, out of the total 
51-ks exposure. In addition, Optical Monitor (OM) images 
were taken in the UVW1, UVM2 and UVW2 filters, with 
exposure times of 3.0, 6.0 and 17.8 ks, respectively 2 . We fil- 
tered the EPIC event files, selecting only the best-calibrated 
events (pattern ^ 12 for the MOSs, pattern ^ 4 for pn), and 
rejecting flagged events. 

We used the Chandra Interactive Analysis of Observa- 
tions (CIAO) task wavdetect to identify discrete sources in 
the combined EPIC image, in different energy bands 3 . We 
selected a source extraction region of radius 25" around the 
brightest point-like source (the only one of the two ROSAT 
ULXs that is seen again by XMM-Newton, see Section 2.2). 
We built response functions with the SAS task rmf gen, and 
auxiliary response functions with arf gen, for pn and the two 

2 See |http:// www.xmm.ac. uk/onlines /uhb/XMM_UHB/node66 . html 
for a plot of the OM filter throughput curves. 

3 A comparison between the source-finding algorithms in CIAO 
and SAS can be found in Valtchanov, Pierre & Gastaud (2001). 



MOSs. We then coadded the three spectra, with the method 
described in Page, Davis & Salvi (2003), in order to in- 
crease the signal-to-noise ratio. Similarly, we obtained com- 
bined EPIC spectra of the X-ray emission from the whole 
galaxy (excluding the ULX) and from the star-forming re- 
gion within V from the galactic centre. Finally, we fitted 
the background-subtracted spectra with standard models in 
XSPEC version 11.3.1 (Arnaud 1996). 

2.2 Discrete X-ray sources and unresolved 
emission 

The X-ray emission from M 99 is dominated by hot g 
expected from its high star formation rate (SFR) . The mod- 
est spatial resolution of XMM-Newton and the strength of 
the diffuse component make it difficult to resolve individual 
point-like sources (essentially, high-mass X-ray binaries and 
young supernova remnants), especially in the inner arcmin 
(4.8 kpc at the assumed distance), where most of the X-ray 
emission is concentrated (Figure 1). A few point-like sources 
are instead detected (Table 1) from the combined EPIC im- 
age in the 1.5-12 keV band, where the diffuse emission is 
negligible 4 . The location of the detected sources is shown 
in Figure 2, overplotted on a true-colour OM image. Of the 
two ULXs found by ROSAT inside the D 25 ellipse in 1997, 
one (X-l in the ROSAT catalogue, Liu & Bregman 2005) is 
no longer detected: hence, we estimate that it must now be 
fainter than » 5 X 10 38 erg s" 1 in the 0.3-12 keV band; its 
luminosity was m 6 x 10 39 erg s _1 in 1997, extrapolated to 
the same band. The other ULX (IXO 46 in Colbert & Ptak 
2002, or X-2 in Liu & Bregman 2005) is instead recovered, 
with a 0.3-12 keV luminosity « 2 x 10 40 erg s~ x (detailed 
analysis of this source in Section 2.4). 

We estimate that the total unabsorbed X-ray luminosity 
(not including the bright ULX) from inside the D25 ellipse 
is « (1.2 ± 0.2) x 10 41 erg s -1 in the 0.3-12 keV band, of 
which only » 10 40 erg s" 1 in the 2-12 keV band. About 
70% of the X-ray emission comes from within a projected 
radius of 1' ~ 4.8 kpc. The unresolved emission is dom- 
inated by thermal plasma at a temperature < 0.30 keV; a 
two-temperature thermal-plasma fit (vapec in XSPEC) gives 
temperatures kT x = (50 ± 5) eV and kT 2 = (0.30 ± 0.01) 
keV (Figure 3). Gas temperatures < 0.3 keV are typical 
of disk emission in normal spirals, and are also one of the 
plasma components in starburst galaxies. We do not detect 
any additional thermal-plasma component at ~ 0.7-0.8 keV, 
which is instead characteristic of starbursts (e.g., Ott, Wal- 
ter & Brinks 2005; Jenkins et al. 2005; Summers et al. 2003, 
2004). The diffuse, soft X-ray emission, which traces recent 
star formation, appears to be uniformly distributed over the 
inner disk, and only weakly peaked in the nuclear region. 

An additional, harder component of the unresolved 
emission is well fitted by a power-law of photon index 
r = 1.70 ± 0.25. As is generally the case for star-forming 
galaxies, the power-law component can be interpreted as 
emission from unresolved X-ray binaries. The contribution 
of the discrete sources (both resolved and unresolved, but 
not including the ULX) to the total emitted luminosity in 

4 The same five sources listed in Table 1 are also found when we 
use the SAS source-finding routine emldetect. 
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No. 


R.A.(2000) 


Dec. (2000) 


/l.5-12 
(10 14 ^rg cm 2 s l) 


/o.3-12 
(10 14 ^rg cm 2 s ) 


-^0.3-12 

/ 1 fl39 a-tre o ~ 1 "\ 

erg & ) 


01 


12 18 45.58 


14 26 43.1 


1.71 ±0.35 


2.46 ±0.50 


0.85 ±0.20 


02 


12 18 51.13 


14 26 30.4 


2.09 ±0.38 


3.01 ±0.55 


1.0±0.2 


03 


12 18 51.36 


14 24 24.5 


3.02 ±0.51 


4.35 ±0.73 


1.5 ±0.3 


04 


12 18 52.61 


14 25 47.8 


2.44 ± 0.38 


3.52 ±0.55 


1.2 ±0.2 


05 


12 18 56.15 


14 24 18.0 


39.5t 2 -.2 


56.]+?;* 


19 4+ 1 ' 2 



Table 1. Point-like X-ray sources detected at > 3-cr significance in the 0.3-12 keV band inside the D25 ellipse of M99. Fluxes and 
luminosities are isotropic, emitted values, assuming a power-law spectrum with V = 1.7 and a line-of-sight Galactic column density 
A^.Gai = 2.7 x 10 20 cm" 2 . 




Table 2. Best-fitting parameters to the combined XMM- 
Afeiuion/EPIC spectrum of the inner-disk emission; the source 
region is a circle of radius 1'. The XSPEC model is tbabSQ a i X 
tbabs X (powerlaw + vapeci + vapec2). The quoted errors are 
the 90% confidence limit and A^Gal = 2.7x 10 20 cm~ 2 (Dickey & 
Lockman 1990). K\ and K 2 are the normalisation constants of the 
two vapec components; fS h %_ 12 is the observed flux; /q™_ 12 the 
unabsorbed flux; (L po / L)o.3-i2 the relative contribution of the 
point-like accreting sources (resolved and unresolved) to the total 
X-ray emission, clearly dominated here by hot thermal plasma. 



Parameter Value 



N H (10 20 cm" 2 ) 




r 




K po (10- 5 ) 


0+1.0 

°-°-1.2 


fcTi (keV) 


o.05ol° ;» 


K x (lO" 3 ) 


3501??° 


kT 2 (keV) 


S0+ 01 


K 2 (lO" 3 ) 


1 5+ 1 ' 1 
1 -°-0.2 


(Z/Zq) 


13+ 036 


xl 


0.96(81.0/84) 


/0.3-12 (10- 12 CGS) 


30+0 02 


foS-a ( 10 " 12 CGS ) 


n r+2.4 
Z -°-1.5 


io.3-12 (10 40 ergs- 1 ) 




(Lpo/i)o.3-12 


« 0.10 



the 0.3-12 keV band is only ~ 15% of the total luminos- 
ity inside the D25 ellipse. The relative importance of the 
thermal-plasma emission over the point-source contribution 
is even higher in the inner-disk region (Table 2), where the 
power-law component is only « 10% of the total luminosity. 



2.3 Nuclear morphology: an old conundrum 

The nucleus itself is not unequivocally resolved as a point 
source. We do clearly detect X-ray emission in the hard band 
at the nuclear position (Figure 1, middle panel), for an es- 
timated emitted luminosity w 10 39 erg s _1 extrapolated to 
the 0.3-12 keV band. However, because of the low number 
of counts and low spatial resolution, we cannot tell whether 
it is the true nuclear source, or instead a number of unre- 




Figure 2. True-colour UV image from XMM-Newton/OM.; red 
corresponds to the UVWl filter, green to the UVM2 filter, and 
blue to the UVW2 filter. Point-like X-ray sources detected with 
XMM-NewtonfEPIC are marked with green circles, of radius l'/5. 
The peak of the X-ray emission in the nuclear region is indicated 
with a red circle, of radius 3" . 

solved X-ray binaries in the region where star formation is 
more active. 

On an historical note, the nature of the nuclear source in 
M 99 was the subject of intense investigations by Carl Lam- 
pland at the Lowell Observatory, many decades ago (Lam- 
pland 1921). Lampland noted that the peak of the optical 
emission in his photographic plates (which have peak sensi- 
tivity at w 4200A, in the standard B band) appeared some- 
times shifted from the position corresponding (as we know 
today) to the true galactic nucleus, to another point located 
w 2'.'5 to the southeast. This variability was re-examined and 
confirmed by Walker (1967), who found that the off- nuclear 
source was brighter in 10 of the 30 Lowell plates taken be- 
tween 1921 and 1948, with changes occurring sometimes on a 
timescale of a few days. Spectroscopic observations (Walker 
1967) also suggested that the true optical nucleus is redder, 
corresponding to a stellar type G0-G2, while the off-nuclear 
source appears much bluer. 

We can now shed further light on the nuclear morphol- 
ogy. An HST/WFPC2 image in the F606W filter (Figure 
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Figure 1. XMM-Newton/'EPIC X-ray maps of the inner part of M99. Left panel: combined EPIC image and contours for the soft band 
(0.2-1.5 keV), dominated by emission from diffuse hot plasma. Middle panel: EPIC image in the hard band (1.5—12 keV), with soft-X-ray 
contours superimposed for comparison. Right panel: greyscale optical image in the Bessell- il band (archival data from the VXT/FORS1) 
with soft-X-ray contours overplotted. 



4, left panel) is dominated by a bright nuclear star cluster 
(a nuclear structure typical of late- type spiral galaxies), lo- 
cated at the dynamical centre, as can be inferred from the 
morphology of the dust filaments. The cluster is spatially 
resolved in the HST image, with a full width half maxi- 
mum of 0'.'15 ± 0'.'02 corresponding to (12 ± 2) pc; point 
sources in the same field have a width of O'.'lO ± O'.'Ol. An 
approximate conversion to standard colours gives an ab- 
solute magnitude My ps —13.5 mag. Within 200 pc from 
the nucleus, there are a few, fainter, unresolved stellar clus- 
ters with typical absolute magnitudes My w — 11 mag. The 
secondary optical peak described by Lampland (1921) and 
Walker (1967) corresponds to the stellar complex labelled 
as A in Figure 4 (unresolved from ground-based telescopes) , 
and sometimes (at least in plates with worse seeing) to the 
unresolved blend of A and B. When we compare the HST 
image with an XMM-Newton/OM image of the same field, 
in the near-UV (UVW2 filter, with an effective wavelength 
of 2120A), we note that the true nucleus is undetected: the 
brightest UV sources in the nuclear region are the young 
stellar complexes labelled as A, B and C (Figure 4, right 
panel). The optical nucleus remains undetected even in the 
UVW1 filter (effective wavelength of 2910A). Unfortunately, 
the XMM-Newton/EPIC image does not provide enough in- 
formation to determine whether the X-ray emission peaks 
at the true optical nucleus or at any of those three UV- 
bright sources nearby. It is plausible that the true nucleus 
and source A would appear of comparable brightness in the 
standard U band. The nuclear cluster is the brighter source 
in the B band, as evident from the archival images in the 
NASA/IPAC Extragalactic Database 5 (NED). Of course, 
this does not explain how the peak brightness could be per- 
ceived to shift between the two sources several times over 
a few decades, in plates with the same density, taken a few 
days apart. Repeated supernovae or variability of individual 
stars are clearly unsatisfactory explanations; variations in 
the brightness of a nuclear BH would also be very peculiar. 
The true explanation of this effect, if real, remains a mystery 
for now. 



5 http:/ /nedwww. ipac.caltech.edu/ 
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Figure 3. Coadded XMM-Newton/EPIC spectrum and best-fit 
residuals of the unresolved emission in the inner disk (within 1' 
from the nucleus). The spectrum has been fitted with a two- 
tempcraturc thermal plasma component plus a power law. See 
Table 2 for the best-fitting parameters. 



2.4 A bright, hard ULX 

The bright ULX located w 8 kpc southeast of the nucleus 
was the main objective of our X-ray study. Above 1.5 keV, 
the ULX is brighter than the rest of the galaxy (Figure 1). Its 
background-subtracted X-ray spectrum (~ 2100 net EPIC 
counts) is well fitted by a simple absorbed power-law (Fig- 
ure 5 and Table 3) of index T — 1.7 ± 0.1. There is no hint 
of a thermal component (blackbody or disk-blackbody) that 
could be attributed to an accretion disk. Its emitted lumi- 
nosity in the 0.3-12 keV band is l-9±ol? x 10 40 erg s" 1 , which 
is the same, within the errors, as the extrapolated ROSAT 
luminosity observed in 1997 (Liu & Bregman 2005). Inter- 
estingly, the X-ray emission detected by Einstetn/lPC (Fab- 
biano, Kim & Trinchieri 1992; archival X-ray image avail- 
able in NED) is centred on the galactic disc and does not 
show an enhancement at the ULX position, even allowing 
for the lower spatial resolution and low number of counts 
(140 ± 16). As a further comparison, there is another X-ray 
source clearly visible in the Einstein image, about 8' west 
of the galaxy (a background quasar). In the XMM-Newton 
observation, its flux in the Einstein energy band is less than 
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Figure 4. Morphology of the nuclear region in the optical and near-UV. Left panel: ffi?T/WFPC2 image in the F606W filter; right 
panel: the same region, seen by XMM-Newton/OM in the UVW2 filter. At the dynamical center, an old nuclear star cluster dominates 
the optical image but is undetected in the near UV. The labels identify smaller, younger clusters of OB stars, and are given simply to 
facilitate a comparison between the two images. 



half of the ULX flux; however, that background source was 
clearly stronger than the ULX in the Einstein image. We 
estimate that at the time of the Einstein observations (1980 
June 25), the ULX must have been at least a factor of 5 
fainter than in 1997 and 2003. 

We estimate that the XMM-Newton/EPIC astrometry 
is accurate to ~ l'.'5. There are no other point-like X-ray 
sources in the field that can allow a better registration 
with optical positions. We searched for optical counterparts 
within the X-ray error circle in the XMM-Newton/OM im- 
ages and in several optical images from public archives but 
found none. The bright optical source suggested as a pos- 
sible counterpart by Colbert & Ptak (2002), based on the 
ROSAT/HUl position, is now clearly ruled out (Figure 1, 
right panel). Its brightness and colours are consistent with 
a GO foreground star (Kharchenko 2001). We took optical 
spectra of this source with the LRIS spectrograph on the 
Keck telescope, confirming the spectral identification and 
verifying that its radial velocity is consistent with a fore- 
ground star, ruling out the possibility that it belongs to 
M 99. Using the best high-resolution optical images available 
in public archives (in particular, a 30-s VLT/FORSl obser- 
vation in the Bessell-i? band; a 720-s B-band image from the 
2.1-m telescope at Kitt Peak National Observatory; a 200-s 
B-band image from the 2.5-m Isaac Newton Telescope at La 
Palma), we conclude that there are no optical counterparts 
brighter than Mb ~ —9.5 mag and Mr « —9 mag. This is 
not a very stringent limit, and does not even rule out an old 
globular cluster. None the less, we can confidently say that 
the ULX is right at the outer edge of the stellar disk (the 
limit beyond which the gas density is too low to collapse 
and form stars spontaneously), and that there are no super 
star-clusters or massive OB associations at that position or 
within ~ 800 pc. 



Table 3. Best-fitting parameters to the combined XMM- 
Newton/EPlC spectrum of the ULX. The XSPEC model is 
tbabSQ a i X tbabs X po. The quoted errors are the 90% confi- 



dence limit and iVjj Gal 
1990). 



2.7 X 10 20 cm _i (Dickey & Lockman 



Parameter 



Value 



N H (10 20 cm" 2 ) 

r 

K po (10- 5 ) 



23.3 
1.67 
7.3 



+3.8 
3.4 
+ 0.10 
0.10 
+0.9 
0.8 



obs 
3- 



f ol 
JO. 

/em 
0.3- 

^0.3 



(10~ 13 CGS) 
(10- 13 CGS) 
(10 40 erg s" 1 ) 



0.89(86.0/97) 
4.4 



+0.2 
0.3 



O.D_ 2 



1.9 



+0.2 
0.1 



3 RADIO STUDY: A COLLIDING CLOUD? 

If the optical data offer no clue on the nature of the ULX, 
the radio data may instead suggest a complex, intriguing 
story. H I 21-cm line observations of M 99 were carried out by 
Phookun et al. (1993), using the C and D arrays (exposure 
times of 8 hr and 4 hr, respectively) of the Very Large Array 6 
(VLA). The size of the synthesized beam was 25'.'03 x 23'.'56. 
See Phookun et al. (1993) for details of the observations and 
data analysis. As expected, the H I disk is ~ 40% larger than 
the stellar disk, and extends slightly beyond the location 
of the ULX (Figure 6). The total HI flux from the galaxy 
corresponds to a gas mass ss 5 x 1O 9 M0 (Vollmer et al. 2005; 



6 The VLA is a facility of the National Radio Astronomy Obser- 
vatory, which is operated by Associated Universities, Inc., under 
cooperative agreement with the National Science Foundation. 
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channel energy (keV) 

Figure 5. Coadded XMM-Newton/EPIC spectrum and best-fit 
residuals of the ULX spectrum, fitted with a power law of index 
r = 1.7 ± 0.1. See Table 3 for the best-fitting parameters. 



Phookun et al. 1993), using a standard conversion between 
HI flux and mass (e.g., Rohlfs & Wilson 2006). 

Interestingly, the ULX is located close to where a large 
HI spur or cloud apparently joins the gas disk (Figure 6). 
The gas mass of this cloud is ~ 10 7 Mq, considering only the 
fraction that visibly "sticks out" from the disk. More likely, 
its mass could a factor of two higher, considering that part 
of it overlaps with the disk (see also Fig. 5 in Phookun et 
al. 1993). This cloud represents ~ 10% of the gas mass out- 
side the galactic H I disk: either infalling onto the disk, or in 
the process of being shredded from it, according to alterna- 
tive models (recall Section 1). It may be possible that the 
association is simply a coincidence: there is no definitive way 
to tell from the data available. It is none the less remarkable 
that, while there are various other gas clouds around the 
galaxy (references in Section 1), the cloud near the ULX is 
the only one that overlaps and probably impacts the stellar 
disk; all the others are at larger radial distances, at the edge 
of the galactic H I disk. We estimate that only « 3% of the 
projected area of the stellar disk overlaps with this (or any 
other) external cloud. 

The H I velocity map (Phookum et al. 1993) tells a more 
dramatic picture (Figure 7). The gas cloud is strongly blue- 
shifted with respect to the galactic disk, with a difference in 
the projected radial velocity > 100 km s , and is clearly 
overlapping with the disk. From the inferred radial veloci- 
ties, the cloud could be either in the process of being ejected 
(ram-pressure stripped) from the disk towards us, or, more 
likely in our opinion, it is tidal debris infalling onto the disk 
from behind. In the latter case, we have no elements to de- 
termine whether the cloud is actually hitting and merging 
with the disk. But is is interesting to note that the possible 
impact point is very close to the projected ULX position. 
We shall discuss in Section 4 whether there may be a di- 
rect physical connection between collisional processes and 
the formation of ULXs. 




Figure 6. Top panel: H I 21-cm radio flux, represented in 
greyscale over the galactic disk, and as contours for the outer 
(less dense) regions. Contour levels are on a square-root scale 
from 0.02 to 0.3 Jy km s beam . Bottom panel: the same 
HI radio flux contours overplotted over a VXT/FORSl ij-band 
image. 



4 DISCUSSION 

4.1 Hard state of ULXs 

With an X-ray luminosity « 2 x 10 40 erg s _1 , the ULX 
in M99 is among the brightest of this mysterious class 
of sources: the cut-off in the luminosity function is at « 
3x 10 40 erg s" 1 (Gilfanov, Grimm & Sunyaev 2004; Swartz et 
al. 2004) . If the emission is isotropic and Eddington-limited, 
the mass of the accreting BH is > 100A/©. Apart from its 
extreme brightness, there is a notable, unexplained feature 
associated with this source and other bright ULXs: the X- 
ray spectrum is well fitted by a simple power-law with a 
rather hard photon index (r = 1.7 ± 0.1). In stellar-mass 
BH binaries, such spectra are characteristic of the low/hard 
state (e.g., McClintock & Remillard 2006), with typical X- 
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Figure 7. Velocity map of the HI 21-cm radio emission. Char- 
acteristic radial velocities (in km s" 1 ) are also marked on the 
plot, at various places on the disk and the cloud. The gas cloud is 
consistent with tidal debris infalling onto the galactic disk from 
behind (with respect to our viewpoint), at a projected relative 
speed > 100 km s _1 ; it appears as though it is impacting or 
merging with the disk at or near the ULX position. 



ray luminosities in the standard XMM-Newton or Chandra 
bands less than a few per cent of the Eddington luminosity. 
At higher mass accretion rates and luminosities, stellar-mass 
BH binaries are dominated by either an accretion disk com- 
ponent (thermal, or high/soft state) or by a steep power-law 
component (r > 2.5). If ULXs follow the same spectral-state 
classification as stellar-mass systems, the inferred spectral 
slope and luminosity of the source in M99 would suggest 
an accreting BH with a mass ~ 1O 4 M0. This would require 
so-far untested formation processes, such as primordial rem- 
nants or nuclei of accreted dwarf galaxies that have recently 
captured a donor star. Alternatively, ULXs may simply have 
a different spectral behaviour; for example, they might be 
found in a steady high/hard state which is not known in 
Galactic X-ray binaries. For other examples of bright ULXs 
with relatively hard power-law-like spectra, see Winter et 
al. (2005), and Stobbart, Roberts & Wilms (2006); see also 
the discussion in Gongalves & Soria (2006) and references 
therein. 

The underlying problem is how to explain the lack of 
evidence or the relative weakness of the accretion disk com- 
ponent in such bright sources. If the disk is faint because it 
is truncated (as in the standard low/hard state), we expect a 
low radiative efficiency, as well as a dimensionless accretion 
rate much below Eddington; therefore, a very high accre- 
tion rate (in physical units) and an even higher BH mass 
are needed to explain the observed luminosity. If the disk 
extends to the innermost stable circular orbit, the efficiency 
may be higher (standard efficiency ~ 0.1); however, most of 
the detected photons must have been reprocessed in a comp- 
tonizing medium, because we do not see the thermal disk 
spectrum directly. This requires that most of the gravita- 
tional energy is released in or transferred to the upscattering 
medium, for example via magnetic coupling (Kuncic & Bick- 



nell 2004). Finally, especially for a pure power-law source, 
we cannot rule out that the X-ray emission comes from a 
relativistic jet pointing towards our line of sight (microb- 
lazar) (Kording, Falcke & Markoff 2002). However, physical 
and statistical arguments make the strong-beaming scenario 
rather implausible as a general explanation for the whole 
population of ULXs (Davis & Mushotzky 2004). Direct evi- 
dence against beamed emission is provided in some cases by 
the energy requirements of a photoionized nebula surround- 
ing a ULX (Pakull & Mirioni 2003) . Relativistic beaming has 
also been ruled out, in one case, by the presence of quasi- 
periodic oscillations (M82 X-l: Strohmayer & Mushotzky 
2003), and in another case by the radio/X-ray flux ratio 
(Holmberg IX X-l: Miller, Fabian & Miller 2004). 

4.2 ULXs and galaxy collisions 

Bright ULXs are preferentially associated with tidally inter- 
acting or collisional systems, or actively star-forming galax- 
ies, or both (e.g., Swartz et al. 2004 for a population study). 
For example, many bright ULXs have been found in the 
Antennae (Zezas & Fabbiano 2002), the Cartwheel (Gao 
et al. 2003), the Mice (Read 2003), NGC 7714/15 (Smith, 
Struck & Nowak 2005), NGC 4485/90 (Roberts et al. 2002) 
the M81/M82 group (in M82, Holmberg II and Holmberg 
IX). 

The simplest explanation is that galaxy collisions trig- 
ger or enhance star formation, which, in turn, increases the 
birth rate of high-mass X-ray binaries (hence, the normaliza- 
tion of the point-source X-ray luminosity function). This also 
increases the probability of forming very luminous sources, 
near the upper cut-off. For example, it explains the large 
number of X-ray sources with luminosities > 10 39 erg s _1 in 
the Antennae (Zezas & Fabbiano 2002). 

While this is probably a correct argument, it cannot 
be the whole story, especially for the small sample of ULXs 
with luminosities > 10 40 erg s _1 . Some of them (e.g., Holm- 
berg H X-l: Dewangan et al. 2004; Holmberg IX X-l: Miller 
et al. 2004) are in dwarf galaxies with relatively small star- 
formation rates. Others, such as the ULX in M 99 or the 
brightest ULXs in NGC 7714 (Soria & Motch 2004) and 
NGC 4559 (Soria et al. 2005), are in strongly star-forming 
galaxies but far away from the main star- forming regions. 
There are no ULXs in the inner disk of M 99 despite an 
SFR ~ 10M Q yr _1 ; instead, the bright ULX object of our 
study is at the outer edge of the stellar disk, where the star 
formation rate is and probably has always been orders of 
magnitude lower. And yet, the ULX does seem to be as- 
sociated with a collisional event between the disk and an 
infalling gas cloud (Section 3). How can we make sense of 
this apparently contrasting evidence? Here, we try to specu- 
late one possible scenario for ULX formation, consistent with 
the sketchy observational evidence available so far. The fol- 
lowing statements are conjectures to be tested with further 
observational and theoretical modelling. 

i) Most ULXs are not strongly beamed. They are powered 
by BHs with masses ~ 2O-2OOM0 ; the upper limit can easily 
be reduced to < lOOM© if mild anisotropy or mild super- 
Eddington emission are allowed. The donor is likely to be 
a Roche-lobe-filling OB donor star, which can supply the 
required amount of accreting gas over its nuclear timescale 
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(Rappaport, Podsiadlowski & Pfahl 2005; Copperwheat et 
al. 2006). If so, ULXs represent the upper end of high-mass 
X-ray binaries, consistent with their preferential location in 
young stellar environments; the accreting BHs are a factor 
of a few more massive than Galactic systems, and require 
stellar progenitors with masses ~ 5O-4OOM0. 

ii) A starburst or high SFR may favour but is not a suffi- 
cient condition for the formation of bright ULXs; other con- 
ditions may be required. For example, low metal abundance 
is essential to allow the evolution of a very massive star into 
a massive BH (Pakull & Mirioni 2003; Heger et al. 2003). 
A high SFR is not a necessary condition, either, as proved 
by the occasional presence of ULXs in gas-rich but low-SFR 
environments. Super star-clusters are also not a necessary 
condition for ULX formation. 

in) Galaxy-galaxy or cloud-galaxy collisions may induce 
the formation of bright ULXs, as well as trigger intense star 
formation, as separate, parallel consequences. This could ex- 
plain why the two phenomena (ULXs and starbursts) are 
often, but not always, associated with each other. For exam- 
ple, local collisional events in NGC 4559 (Soria et al. 2005) 
and M99 may be responsible for their respective ULXs in 
metal-poor regions at the very edge of their stellar disks. 

The key issue we need to explain is why a collisional 
event would directly (i.e., not simply through enhanced star 
formation) favour the formation of a very massive star, 
which may then become a ULX progenitor if other condi- 
tions (e.g., low metal abundance) are also satisfied. A possi- 
ble explanation is suggested by studies of the star-formation 
process in the Galactic protocluster NGC 2264-C, in the 
Cone Nebula (Peretto, Andre & Belloche 2006). In that case, 
a molecular clump with a mass « 1700 A/q is undergoing 
a large-scale collapse onto its central region, on a dynam- 
ical timescale (essentially in free fall). This results in the 
formation of a few massive protostars, one of which has al- 
ready reached a mass ~ 40Mq but is still accreting at a rate 
~ 10~ 3 M© yr" 1 (Peretto et al. 2006). It is plausible that two 
or three protostars would accrete a few lOOM© of gas and co- 
alesce even before ending their protostellar phase (~ 3 x 10 5 
yr). Based on this scenario, we speculate that the formation 
of a star with an initial mass of a few lOOM© may occur in 
a protocluster such as NGC 2264-C, via massive global gas 
infall and mergers already during the protostellar phase (So- 
ria 2006); this is perhaps a more common process than the 
runaway merger of O stars in the core of super star-clusters 
(Portegies Zwart & McMillan 2002; Freitag, Giirkan & Rasio 
2006). 

It can be shown that an external shock or pressure wave 
can trigger the global, dynamical collapse of a molecular 
clump such as NGC 2264-C; instead, an isolated clump in 
hydrostatic equilibrium would have fragmented into much 
smaller Jeans-mass protostars (for recent reviews, see for 
example Struck 2005, Elmegreen 2002, 2004). The critical 
momentum required to trigger the collapse of a molecular 
clump is of order of the clump mass times its sound speed; 
moreover, the accretion rate onto the protostellar cores is 
proportional to the momentum imparted to the collapsing 
molecular clump (Motoyama & Yoshida 2003). Supernova 
shocks are a possible trigger of dynamical collapses in nearby 
molecular clumps; a fast, colliding gas cloud may carry even 
more energy and momentum than a supernova. Thus, we 



speculate that collisional events (such as the cloud-disk col- 
lision at the outer edge of M 99, or the possible dwarf galaxy- 
disk collision in NGC 4559: Soria et al. 2005) may lead to 
the formation of extremely massive protostars, via dynam- 
ical collapse and protostellar mergers on a timescale ~ 10 5 

yr- 

The second step of the process to form a ULX would 
then be to ensure that the massive stellar progenitor re- 
tains most of its gas until core collapse. We speculate that 
this is where low metal abundance comes into play, because 
it strongly reduces the mass loss rate through line-driven 
stellar winds (e.g., Eldridge & Vink 2006, and references 
therein) and therefore leaves behind a bigger BH. This is 
why the (metal-rich) Pistol star in the Quintuplet Cluster 
(initial mass « 200-250M Q : Figer et al. 1998), the Wolf- 
Rayet stars in the Arches cluster, or the massive protostars 
in the NGC 2264-C protocluster will never form massive 
BHs, and perhaps why there are no ULXs in that part of 
our Galaxy. Conversely, we expect metal abundance to be 
much lower at the outer edge of the disk in M99, perhaps 
made even lower by the infalling gas clouds (see, e.g., Chen, 
Hou & Wang 2003, and Andrievsky et al. 2003, for recent 
studies of metal abundances as a function of galactocentric 
distance in the Milky Way; and MacArthur et al. 2004 for 
similar trends in other disk galaxies). 

In summary, our suggested scenario differs from models 
based on the runaway coalescence of O stars in super star- 
clusters, because it does not require such massive systems. 
We speculate that the formation of a massive stellar pro- 
genitor may occur in smaller (~ 1O 3 -1O 4 M0) protoclusters, 
during their embedded phase, through infall and merger pro- 
cesses similar to what is currently observed in NGC 2264-C. 
Protoclusters of that size are massive enough to produce O 
stars, but would by no means be classified as super clusters; 
in fact, they may not even be massive enough to survive 
their embedded phase as bound systems (Kroupa & Boily 
2002). 

Our scenario also differs from ULX formation models 
based on Population-Ill stars, because it only requires low 
but not primordial abundances, and therefore it can work 
at any redshift. At zero metallicity, stars with initial masses 
« 14O-3OOM0 may be disrupted by pair-instability super- 
novae that leave no remnants (Heger et al. 2003). Instead, 
at Z ~ O.IZq, stars in that mass range may become the 
progenitors of accreting BHs in ULXs, via direct collapse, 
precisely in the mass range required by the X-ray observa- 
tions. Further discussion of each of those speculations, and 
detailed comparisons with the observations, is beyond the 
scope of this paper. 



5 CONCLUSIONS 

We studied the X-ray properties of M99, the most massive 
spiral galaxy in the Virgo Cluster. As expected from its high 
SFR (« 1OM yr" 1 ), the X -ray emission is dominated by 
a soft thermal-plasma component. The total unabsorbed lu- 
minosity (not including a bright ULX) inside the D25 ellipse 
is w (1.2 ±0.2) x 10 41 erg s _1 in the 0.3-12 keV band; about 
15% of this is due to resolved or unresolved discrete sources. 
The emission appears almost uniformly diffused across the 
inner disk (at < 5 kpc from the nucleus), and there is no 
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starburst core. The temperature of the hot gas, kT < 0.30 
keV, is also rather low, typical of disk emission rather than 
of a starburst environment. 

We briefly discussed the morphology of the nucleus in 
the optical and UV bands, showing the presence of a (redder) 
massive nuclear star cluster and a few smaller, much bluer 
clusters of young stars around it. The spatial resolution of 
XMM-Newton does not reveal whether there are faint point- 
like X-ray sources in the nuclear region, and if so, whether 
they are located in the old nuclear cluster or associated to 
the young stars around it. 

The main goal of this paper was a study of the prop- 
erties and origin of a bright ULX at the outer edge of the 
stellar disk. Its unabsorbed luminosity of ~ 2 x 10 40 erg 
s _1 in the 0.3-12 keV band, together with a pure power- 
law spectrum with photon index F « 1.7, are difficult to 
interpret in the framework of classical spectral states for 
stellar-mass BHs. Such high/hard states are not uncommon 
in ULXs but are not generally seen in Galactic BHs. It sug- 
gests that we are not seeing the accretion disk directly, and 
that most of the gravitational power is efficiently transferred 
to a comptonizing region. 

An intriguing new discovery is that there is a massive 
gas cloud (HI mass ~ 1O 7 M0) seen in projection very close 
to the ULX position. From its radial velocity, we suggested 
that this cloud is falling onto the galactic disk (from be- 
hind) and perhaps impacting it near the ULX location. The 
cloud occupies « 3% of the projected area of the stellar disk; 
therefore, we cannot entirely rule out the possibility that it 
is a chance association. However, we explored the specula- 
tive idea that there is a direct connection between collisional 
events and ULX formation, and not simply an indirect effect 
due to an enhancement in star formation. We suggested a 
possible scenario of ULX formation that would be consistent 
with this interpretation. We argued that cloud collisions may 
trigger a large-scale dynamical collapse of molecular clumps, 
rather than Jeans-mass fragmentation; in turn, this may lead 
to the formation of very massive stellar progenitors in the 
cluster core. If the metal abundance is low (Z ~ O.IZq, not 
primordial abundances), we speculate that the star might 
retain enough of its gas to directly-collapse into a BH mas- 
sive enough (up to ~ IOOMq) to explain the luminosity of 
the brightest ULXs. 

This scenario has to be tested with further individual 
and population studies of ULXs. At the individual level, the 
suggestion that ULXs have masses up to an order of magni- 
tude higher than Galactic BHs (but not more) needs to be 
tested with more advanced X-ray spectral modelling, and a 
better understanding of the thermal disk and comptonized 
emission components in such systems. At the population 
level, we need a more systematic study of the spatial as- 
sociation between ULXs, low-metallicity environments and 
collisional events. Further work on this issue is currenly un- 
der way (D. Swartz et al., in preparation). 
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